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We study the phase diagram of the three flavor Polyakov-Nambu-Jona Lasinio (PNJL) model and 

in particular the interplay between chiral symmetry restoration and deconfinement crossover. We 

compute chiral condensates, quark densities and the Polyakov loop at several values of temperature 

and chemical potential. Moreover we investigate on the role of the Polyakov loop dynamics in the 

1^^^^ ■ transition from nuclear matter to quark matter. 

o 
o 

(N 

> [ I. INTRODUCTION 

o 

^H . A long-standing problem with Quantum Chromodynamics (QCD) is the difficulty of an analytical study in the 

nonperturbative regime. To deal with this issue a few approximate approaches have been developed, based on the 
CN . idea of effective field theories. Among them, the Nambu-Jona Lasinio model [l] (NJL in the following) has become a 

■ popular tool to describe some of the salient aspects of low energy QCD, in particular chiral symmetry breaking and 

(— I I its restoration at high density and/or temperature (see Refs. [4I for reviews). 

O ■ The NJL model neglects the gluon dynamics, the interactions among quarks being described by contact terms. Due 

I I to the absence of gluons, one of the main properties of low-energy QCD at small temperature and baryonic density, 

TT*: namely confinement, is missing in this effective description. This defect not-withstanding, the NJL model enjoys 

(-H ' a significant popularity since it allows analytical treatments in different contexts. In order to cure its deficiencies 

I— I , some extensions have been also proposed. They either try to describe quarks and hadrons in unified approaches (see 

Ref. [S'l and references therein), or add a bag constant to the equation of state of the NJL quark matter [4!]. In the 

latter procedure the deconfinement transition is obtained by computing the pressure of nuclear matter at low density 

and temperature by some effective model and comparing it to the pressures of the quark matter. When the nuclear 

?^ ' matter pressure becomes smaller than the quark pressure, deconfinement is energetically favored. We stress that the 

(T^ r introduction of a bag constant is needed to reproduce the phase transition from the nuclear to the quark phase in this 

Cn ' approach. 

' , Another extension of the NJL model has been suggested in Ref. [5j . In this paper part of the gluon dynamics is 

y—( ' described by a background temporal gluon field coupled to quarks by the QCD covariant derivative. The background 

f — ■ , field adds a potential term l/({^) to the lagrangian. Its value is determined for any temperature T and quark chemical 

^D ■ potential fi by the minimization of the effective potential and depends on the traced Polyakov loop $ [6|] , related to 

^ I the background gauge field. The resulting model is known as Polyakov-Nambu-Jona Lasinio (PNJL) model. 

• ^ ■ As is well known the Polyakov loop in a pure gauge theory is an order parameter of the deconfinement transition [6| . 

p% ' This peculiarity is related to the existence of a discrete symmetry Z3 of the pure gauge action, which is spontaneously 

3 . broken when deconfinement sets in. The Polyakov loop vanishes in the disordered low temperature phase and is 

■ - - ' different from zero in the high T phase. When dynamical quarks are added, Z^ symmetry is explicitly broken and 

$ can no longer be considered as an order parameter. However, as shown by lattice simulations, its behavior as a 

function of T (increasing from zero to non vanishing values when T increases) can still describe the deconfinement 

crossover. One might assume that this happens also when the chemical potential fj, is varied, which should show a 

link between deconfinement and chiral restoration [3, Q . Therefore in the PNJL model, one introduces a Polyakov 

loop dynamics in the NJL model trying to reproduce both chiral symmetry brea king and q uark confinement. 

The PNJL model with two flavors has been extensively studied [I, [13, [HI, III [11 [Mill [11. In Ref. [13] an extension 
to the three flavor model has been considered, by studying the effects of the Polyakov loop on the thermodynamics 
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of the Color-Flavor-Locked phase of high density QCD [18J near its second order transition to the normal phase. The 
aim of this paper is to extend this analysis to moderate densities. In particular we will study the phase diagram of 
the three flavor PNJL model and the transition from nuclear to deconfined quark matter. We shall consider massive 



quarks because, differently from previous work [17[ , mass effects cannot be neglected at moderate values of the quark 
chemical potential. Therefore we introduce bare quark masses and also compute selfconsistently their in-medium 
values in the mean field approximation. On the other hand, since we are mainly interested in the interplay of chiral 
symmetry restoration and deconfinemcnt crossover we neglect the possibility of color superconductivity, which should 
be produced at higher densities 19, 20, 21, 2^]. While chiral symmetry restoration is adequately described by the 
PNJL model, to establish the deconfinement transition we have to consider also nuclear matter, described by an 
effective field theory, so that a comparison of free energies can be made. 

The plan of the paper is as follows. In Section II we summarize the main features of the three flavor PNJL model. 
Section III is devoted to the study of the phase diagram of the model. In Section IV we compare PNJL quark matter 
with nuclear matter. Finally, in Section V we draw our conclusions. 

II. THERMODYNAMICS OF THE THREE FLAVOR PNJL MODEL 

We consider the lagrangian 

^ = ^^fi *^m7^ - "^/ + M7o) V'/ + A + £6 , (1) 

/ 

where the sum is over the three flavors / (= 1, 2, 3 for u, d, s). In the above equation we have introduced the coupling 
of the quarks to a background gauge field A^ = gS^oAa^Ta via the covariant derivative D^ = d^ — lA^] irif is the 
current mass (we assume to„ — nid)- The quark chemical potential is denoted by /.t. The NJL four-fermion and 
six-fermion interaction Lagrangians are as follows [2|: 

8 

A = G^[(V;AaV^)'+(#75AaV^)'] , (2) 

Ce = -X[detv;/(1+75)V'/' + detV;/(l-75)^/'] , (3) 

where A^ are the Gell-Mann matrices in flavor space (Aq = ^2/3 1/) and the determinant is in flavor space as well. 
Working in the mean field approximation the self-consistent equations for the constituent quark masses are 

Mf =mf - AGof +2 K cr/+i cr/+2 ; (4) 

here af ~ {f f) denotes the chiral condensate of the flavor /, and we define a^ = a^, cr^ — a^- We also introduce the 
quark mass matrix M — diag[Mu, Md, Mg]. The gap equation at T = and // = is 

3A// [^ p' 
af = ^ / , dp , (5) 

which depends on the ultraviolet cutoff A. The parameters are chosen as in Ref. [20'|, i.e. 

m„,d = 5.5 MeV , m^ ^ 140.7 MeV , GA^ == 1.835 , KK^ = 12.36 , A = 602.3 MeV . (6) 

By these parameters one gets m^ ~ 135 MeV, tuk ^ 498 MeV, m^, ~ 958 MeV, m,, ~ 515 MeV and /^ ~ 92 MeV. 

Once the lagrangian is specified, the thermodynamic potential at temperature T is obtained after integration over 
the fermion fields in the partition function: 



n=U[T,<^M + '2G Y: 4-4ifa„a.a.-T^|^Trlog^:^i^ 



(7) 



f—u,d.s 

Here tUn = 7TT{2n -\- 1) are Matsubara frequencies. The inverse quark propagator is given in momentum space by 

S-^ = 70 (p° + M - iA^) -j-p-M , (8) 



where A4 = iAq. A most significant difference between tfie NJL and the PNJL model is the presence in the thermo- 
dynamic potential of the gluon contribution U{T, <&, $) describing the dynamics of the traced Polyakov loop 
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in absence of dynamical quarks. We will consider in the following two different choices of the potential. First we 
assume a polynomial expansion in $ and $: 



Polynomial model : 



T* ~ 2 
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($3 _,_ $3 
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where 



b2{T) = ao + ai{^ 
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(10) 



(11) 



This corresponds to the choice = for the parameter 6 introduced in [l3|. Numerical values of the coefficients 
are [9| 

ao = 6.75, ai = -1.95, 03 = 2.625, ag = -7.44 , 63 = 0.75, 64 = 7.5. (12) 

The remaining parameter Tq is the deconfincmcnt temperature of static and infinitely heavy quarks [(^ . This parameter 
depends in principle on the number of flavors and on the quark chemical potential. We shall investigate two cases. 
First, we assume a constant value Tq = 270 MeV, which corresponds to the deconfinement temperature of static 
and infinitely heavy quarks as computed in lattice QCD (see for example [10]). The second case we study is the 
/i-dependent Tq suggested in |13| . 



with Tr = 1770 MeV, ao = 0.304 and 



To = r^e"i/("o''(^)) , 



6(M) = ^(llA^c-2Ar,)„i^iV,|! 



(13) 



(14) 



(this dependence is motivated by the use of hard dense loop and hard thermal loop results for the effective charge 
111, Hi). The prescription in ^ gives Tq = 178 MeV at /i = 0. 
The second model we consider uses instead of (jlOl) the following logarithmic form 



Logarithmic model 



Z^(T,$,$)_ hiT) 



rp4 



with 



62 (T) = So +ai 



To 



0.2 



$$ + 6(r) ln[l - 6$$ + 4('I>3 + $^) - ($$)2] 



rri \ ^ / 'T^ \ "J 

, b(T) = 63 



Numerical values of the coefficients are as follows f.lO| 

^0 = 3.51, hi = -2.47, h2 = 15.2, 63 = -1.75 



(15) 



(16) 



(17) 



In this second case we only consider the case of a fixed value of Tq: To=183 MeV, corresponding to the Nf = 3 value 
of the deconfinement temperature estimated in Ref. |13j . 

Evaluation of the trace gives in both cases for the thermodynamic potential the expression 

17 = iY[T,$,$] + 2G Y. <^} ~ ^KauCTdas - 6 ^ j ^Efe{K~\p\) 



f—u,d,s 



f—u.d,s ' 



_2r y //4log[l + 3<I>e-'^(^/-^)+3l>e-2'3(^/-M) 
^ J 27r3 ^L 



-2T 5: 

f—u,d.s 



(27r)3 
(27r)3 



-3f3iEf-,i) 



log 



1 + 3$e-'3(^/+^) + 3$e~2''(^^+^) + e-^/JC^^z+A") 



(18) 



with Ef = Jp^ + M? and Mj given by Eq. (|4]) . By searching the global minimum of 17 one can get quark condensates 
af and Polyakov loop for any values of the parameters /x and T . 



III. NUMERICAL RESULTS 

In Fig. [T] we show the phase diagrams of three cases mentioned above (two cases for the polynomial model, one 
for the logarithmic model). They are obtained as follows. For any value of the quark chemical potential, the critical 
temperature is identified with the temperature at which the up quark chiral condensate (7„ has a jump (first order 
transition) , or the derivative of au with respect to the temperature is maximum (second order transition) . The phase 
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FIG. 1: Phase diagram of the PNJL model. The three lines correspond to the following models: To fixed at the constant 
value 270 MeV and polynomial form for the Polyakov potential W($) (highest); To = 183 MeV and logarithmic form of the 
potential (middle); Tq = To(^) and polynomial form (lowest). A solid line denotes a first order chiral transition, a dashed line 
a cross-over for the condensate < uu > and a bold point a critical endpoint (CEP). 



diagrams of the polynomial model correspond to the highest and lowest lines, while the intermediate line is obtained 
by the logarithmic potential. The two lines of the polynomial model are obtained with the fixed value Tq = 270 
MeV (upper line) and with Tq{^) of Eq. (fT5)) (lower) respectively. Qualitatively the three curves are similar, which 
shows the physical results are insensitive to the details of the model. The critical endpoints (CEPs) separating the 
cross-over from the first-order lines are located, for the three cases, at (/X£;,Te) « (300,140) MeV (highest curve), 
{Pe,Te) « (300,115) MeV (middle curve) {he,Te) « (310,68) MeV (lowest curve). At ^ = we find for the three 
cases the cross-over temperatures T^{n — 0) ^ 223 MeV (highest curve), T-^(/i = 0) « 200 MeV (middle curve), 
T^{^ = 0) w 189 MeV (lowest curve). 

The numerical values of the cross-over temperatures at /i = are in the same range of the results found in QCD-like 
theories [2J, [2^, [2^ or in lattice QCD [23| . On the other hand, the CEPs are located at values of ^ much higher than 
the respective values found in QCD /like theories [2J, [2^[20|, lattice QCD [27|] and empirical analysis of the ratio of 
shear viscosity to entropy density [2^|, which suggest there is a CEP at T « 165 MeV and /i « 50 MeV. Clearly the 
details of the models matter; in relation to this problem, in Ref. [l^ an extension of the PNJL model (in the case of 
two flavors) has been proposed using also a C'(-0V')^ interaction term. It was found that the effect of the new kind of 
interaction is to locate the critical end point at higher temperature and lower chemical potential than the NJL/PNJL 
results. One might expect similar effects in the three flavor models considered in the present paper. 

The phase diagrams in Fig. [1] are based on the results of Fig. [2] that shows the T— dependence (at fixed /i) of the 
chiral condensates. The u chiral condensates are represented by dashed lines, those of the strange quark by dot-dashed 
lines (solid lines represent the Polyakov loop $). We have reported results only for the polynomial model for the two 
cases of fixed Tq — 270 MeV (upper panels) and Tq{ii) (lower panels), since the results for the case of the logarithmic 
potential are qualitatively similar to those of the case of the polynomial potential with To = 270 MeV. On the right 
panel of Fig. [2] we have results at /i = 350 MeV, on the left those at /i = 150 MeV. 

From the diagrams of Fig. [5] we can identify another relevant transition temperature: T^, i.e. the temperature 
corresponding to the inflection point of $ (defined as the temperature where d^/dt has a maximum). As discussed 
in the introduction T$ cannot be immediately interpreted as the deconfinement temperature because Z^ is not a 
symmetry of QCD with dynamical quarks. In fact in the next section we shall adopt a different method to study 
the deconfinement transition, based on the comparison of the pressures of the two competing states of matter (quark 
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FIG. 2: Upper panels: Polyakov loop "l> (solid line) and chiral condensates for u quark (dashed line) and s quark (dot-dashed 
hne) versus temperature in MeV, at /i = 150 MeV (left) and /x = 350 MeV (right), for the PNJL model with fixed To = 270 
MeV. For each value of /.i the condensates are normalized to their values at T = 0. Lower panels: Same quantities for the model 
with To = To(m). 



and nuclear) . Nevertheless T$ should approximately coincide with the deconfinement transition, if the guess on the 
relevance of the rise of Polyakov loop for the deconfinement of dynamical quarks is correct. One can note a different 
behavior between the two cases Tq = 270 MeV (fixed) and To(/x). In the former case T$ numerically coincides either 
with the chiral crossover (at low [i) or with the chiral first order phase transition (high /i). In the latter case this 
equality is lost, the shift between T^ and Ty^ being of the order of 30 MeV for \x in the interval [0,320] MeV, and 
shrinking as [i is increased above 320 MeV. Another peculiarity of the To(yu) case is that at high temperatures $ > 1. 
This is not acceptable since $ is defined as the normalized trace of a S'f7(3) matrix and should therefore lie in the 
interval < |$| < 1. The fact that |$| can be larger than one in the PNJL model with a polynomial potential 
is well known in the two flavor case [9|, and to solve this problem several improvements have been suggested, e.g 
the use of the logarithmic potential [5|, [lOj or the interpretation of the Polyakov loop potential as a random matrix 
model ill, 29, 30]. In our numerical calculations within the logarithmic potential model we have checked that |$| is 
always less than one. 

Another possible signature of the deconfinement transition is offered by the behavior of the quark number densities. 
In Fig. [3] we show the results for the scaled quark number densities fif. For each flavor /, n/ is defined as follows: 



nj_ 



1 dn 



(19) 



We plot fif as a function of T at ^ = 150 MeV (left) and ^ = 350 MeV (right) for the two cases Tq = 270 MeV and 
To^fj,). At fixed quark chemical potential the quark number density is almost vanishing below the chiral transition 
temperature, rising quickly in proximity of the transition itself. This behavior has been interpreted in the two fiavor 
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FIG. 3: Upper panels: scaled number densities of quarks njT^ versus temperature (MeV) at /i = 150 MeV (left) and /i — 350 
MeV (right) calculated in the PNJL model with fixed To = 270 MeV; solid lines correspond to up and down quarks, dashed 
lines to the strange quark. Lower panels: the same quantities for the model with To = To{ii). 

case as a simulation of the confincd^dcconfined transition in the PNJL model Q . 

IV. COMPARISON WITH NUCLEAR MATTER 



We shall try here to confront the previous calculations with a concrete phenomenological model for the transition 
from nuclear matter to quark matter. The model will be based on a rather rough modelisation, but we shall see that 
the emerging pattern will still not be very far from that suggested in our previous analysis. In order to locate more 
accurately the deconfinement transition, as well as to study its nature, we now examine the pressure of a gas made 
by hadrons at small temperature and small baryonic chemical potential /is = 3/i and compare it to the pressure of 
the PNJL model, the transition occurring when the difference between the two pressures vanishes. Nuclear matter 
will be described by an effective model based on the non-linear extension of the original Walecka model [31| due to 
Boguta and Bodmer [32, l33l| (WBB in the following) . In its two-flavor version this model allows to predict both the 
compression modulus and the effective nucleon mass at saturation in agreement with their empirical values. The 
spinor content of the model we consider here is the whole lightest baryon octet. We will work assuming SU{3) flavor 
symmetry for the couplings, but including mass differences in the baryon octet. 

The effective Lagrangian of nuclear matter in the WBB model describes the baryons coupled to the a scalar meson 
and the uj vector meson, and takes into account cubic and quartic self-interaction terms for the scalar field: 



^ = E^- 



B [ij'"idf, + igujto^) - (Mb - g^cr)] i^B + -z {d'^adf.a - mla'^) 



+ 2^^^^^^^ - ■^'^t.i'^'"" - 3&M(.9^cr)3 + -cig^a)'^ 



(20) 



In this equation W/^^ 



d^OJiy 



diyUjfj, is the field strength tensor of uj. The fields for the a and w-mesons are denoted 



by cr, cjp respectively and ipB is the Dirac spinor for the baryon B having mass Mb (the sum over the index B is 
over the entire octet). Moreover g^, g^^ are the dimensionless coupling constants of the a and lo fields; 6m and c are 
parameters; Too-, to.^^ are the masses of the cr, w-mesons respectively. Values of the other parameters are in Table [J 





bivi/MN 


c 


iga/m.f (W) 


{g^/m^f (W) 


Set 1 


1.46 X 10"^ 


-1.24 X 10"^ 


9.9262 


4.233 


Set 2 


2.95 X 10"^ 


-1.07 X 10"^ 


11.79 


7.15 



TABLE I: Two sets of parameters for the nuclear matter effective lagrangian, see Eq. (I20|l . We present results for the dimen- 
sionless constant &m/Mjv, related to the cubic interaction couphng constant 6jv/ by the nucleon mass Mjv. 



We use two parameter sets. Set 1 corresponds to the softest nuclear equation of state discussed in [33|; Set 2 to 
the stifFest one. These parameters are identical to those used in the model with nucleons only (without hyperons). 
As discussed in [33| the use of the same parameters is justified since they are obtained by the fit of some nuclear 
properties at the saturation density and zero temperature; in these conditions the hyperons do not play any role and 
cannot infiuence the result. 

We limit ourselves to the mean field approximation, in which the meson field operators are replaced by their 
expectation values; moreover, the ground state of the system we are interested in consists of static and uniform matter. 
This fact implies that the expectation values of the meson fields do not depend on time and space coordinates. At 
finite temperature the equations of motion for meson and nucleon fields are [33| 



QaO- 






Mb - .9<jcr 



Vfc2 + {Mb ~ gaof 



/(eB(k), T) dk - hM{9aCjf - c{g„cjf 






m^ujk = 0, 



b^.{k^ - g^u;'^) - {Mb - ga<j)] i>B{k) = , 
where p is the total baryon density given by 

1 



(21) 

(22) 
(23) 

(24) 



/(eB(k),T) 



by 



1 + exp 



/ ea(k)-M A 



is the Fermi distribution function, and es (k) is the fermion dispersion relation given 



eB(k) = g^ujQ + ^Jk"^ + {Mb ~ gaCrf 



(25) 



The pressure of nuclear matter is given by 

P 

' -y 



-^^M{ga(^f - ■^c{g„(Tf - i^mla'^ + -mlul 



Stt^ 



fc2 



Vfc2 + {Mb - ga<jy 



f{eB{^),T) edk 



(26) 



For the evaluation of the nuclear matter pressure at a given chemical potential and a given temperature we solve 
self-consistently Eqs. pT|) and (l24l) . obtaining the values of g^-a and p that are used in Eq. ((26| . 

When we compare the pressures of nuclear matter and quark matter, we face a normalization problem. As a matter 
of fact, a.t T = p — the self-consistent solutions of Eqs. (PT|) and ([M)) are p = and gcrcr = 0, which imply wq = 
by virtue of the equation of motion Eq. (f^^ . and the vanishing of the phase space integral in Eq. (pS)) . Clearly 
quark pressure must be lower than the nuclear one at relatively low values of density and/or temperature since in 
this region of parameters one is in the confined regime. The quark pressure, which is obtained by the thermodynamic 
potential Q in Eq. (jlSp changing its sign, does not satisfy such a requirement. This problem is well known in NJL 
models |j|, and the usual procedure is to subtract a positive constant B, the bag constant, from the quark pressure. 
The choice of B is not unique, and several different suggestions have been analyzed in the literature Q. We choose 
to fix the bag constant by imposing that the transition temperature from nuclear to quark matter at zero chemical 
potential coincides with the crossover temperature T,j of the PNJL model discussed in the previous Section, i.e. the 
temperature where d^/dt has a maximum. In Table HIl we show the values of the bag constants and of r$ obtained 
for the various models. For each of the quark matter models considered here, the values of B obtained in comparison 





B'/* (MeV) 


T* (MeV) 


fic stiflt (MeV) 


fic soft (MeV) 


Polynomial, To = 270 MeV 


428.8 


223 


543 


617 


Polynomial, To = ro(/i) 


426.1 


160 


500 


588 


Logarithmic 


425.6 


153 


490 


566 



TABLE IL Bag constants, deconfinement temperatures at /^ = and deconfinement chemical potentials at T = for the 
different quark matter models. Values of B obtained in comparison with soft and stiff nuclear matter differ of some part per 
thousand, therefore we show only one value of B for each quark matter model. 
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FIG. 4: Nuclear matter/quark matter phase diagrams in the PNJL model with To — 270 MeV (left) and To{iJ.) (right). Nuclear 
matter contains hyperons and is described by the stiff equation of motion. Nuclear matter is favored in comparison with quark 
matter in the region below the line. 



with soft and stiff nuclear matter differ of some part per thousand, therefore we show only one value of B for each 
quark matter model. 

The results for the two cases are reported in Fig. 31 It shows the phase diagram obtained by comparing nuclear 
(with hyperons) and quark matter. The latter corresponds to the PNJL model with Tq = 270 MeV (on the left panel) 
and Tq = To(fi) (right panel). We stress that a more realistic treatment would require the introduction of different 
coupling constants in the baryon octet, but this is beyond the scope of the present study. The nuclear matter is 
chosen to be the one with the stiff equation of state (Set 2 of Table [l|. The transitions are of the first order. At 
T = we find the critical quark chemical potential /Xc ~ 543 and « 500 MeV for the two cases. For the case of 
the logarithmic potential we find /j,c ~ 490 MeV. If we replace the stiff nuclear matter with the soft one, the phase 
diagram is qualitatively the same as before; in this case we find fic ~ 617 MeV, Hc ~ 588 MeV, fic ~ 566 MeV, 
respectively for the polynomial potential with Tq = 270 MeV, the polynomial potential with Tq = To(/i) and the 
logarithmic potential. 

A comparison with Fig. [1] shows that at small chemical potential the deconfinement temperature almost coincide 
with the chiral restoration temperature for the model with fixed Tq (left panel in FigH]), while for the /i-dependent 
To (right panel in FigU) the two temperatures differ. Moreover, at small temperatures the critical baryonic chemical 
potential is higher for the deconfinement transition. However the values of the chemical potential found here are such 
that one of the color superconductive phases might be energetically favored 19, 20, 21, 22], a possibility we have not 
considered in the present paper for simplicity. 



CONCLUSIONS 



We have studied the phase diagram of the PNJL model with three flavors, considering different versions of the 
model, i.e. with polynonfial and logarithmic forms of the Polyakov loop potential, in the former case with two forms 
of the reference temperature: Tq fixed (=270 MeV) and ^-dependent To(/i). One result is related to the chiral 
symmetry phase transitions and is plotted in Fig. [1] In general we find at /i = crossover temperatures higher than 
the analogous result of the NJL model, namely T — 175 MeV. Moreover we find that the critical endpoint slightly 
depends on the choice of the model, but for all the cases considered in this paper the value of fiE is higher than the 
result obtained by lattice QCD or QCD-like theories {fiE ^ 50 MeV). The PNJL result might be improved by adding 



a new interaction term 0{'ipip)^ in the quark lagrangian as in Ref. [T4I. We leave this investigation to a future project. 

We have also studied the role of the Polyakov loop in the nuclei— >quarks transition. In order to describe nuclear 
matter we have adopted the improved Walecka model with a self-interacting scalar field a and a vector field w^. We 
have chosen two sets of parameters of the nuclear matter lagrangian, corresponding respectively to a soft and to a 
stiff equation of state i3j||. We find similar results in the two cases. The result of this analysis is summarized in 
Fig. [4] (stiff equation of state). The quark pressure is normalized by imposing that the deconfinement crossover of 
the PNJL model and the nuclei^quarks transitions occur at the same temperature at fi — 0. We get that chiral 
symmetry restoration temperature and deconfinement temperature almost coincide at T = in the case of fixed Tq 
and polynomial form of the Polyakov loop potential, while for the /x-dependent Tq the two temperatures differ. We 
also find that at small temperatures the critical chemical potential for deconfinement is higher than the one for the 
chiral transition, see Figs. [1] and S) 

This work might be improved in several different ways. For example one should impose electric and color neutrality 
in the quark matter sector. Also the possibility of a color superconductive phase at high density should be taken 
into account. Moreover the nuclear matter equation of state should be replaced by a more sophisticated one taking 
into account population imbalances due to /3— equilibrium and electric neutrality. We leave all these issues to future 
projects. We finally note that while preparing this text we became aware of an almost simultaneous and independent 
study of the three flavor PNJL model [3J| that has some overlap with our work. 
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